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(up to a translation) of extremal entropic Rankine-Hugoniot discontinuities for systems of conser- 
vation laws (typically 1-shocks, n-shocks, 1-contact discontinuities and n-contact discontinuities 
of large amplitude) among bounded entropic weak solutions having an additional trace property. 
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1 Introduction 

In this article, we develop a theory for uniqueness and global L 2 stability of extremal entropy- 
admissible Rankine-Hugoniot discontinuities (typically 1-shocks, n-shocks, 1-contact discontinu- 
ities and n-contact discontinuities) for a wide class of systems of conservation laws endowed with 
a convex entropy. The uniqueness and stability is shown in the class of bounded weak entropic 
solutions verifying the following trace property. 

Definition 1. Let U G i°°(R + x M). We say that U verifies the strong trace property if for any 
Lipschitzian curve t — > X(t), there exists two bounded functions U-, U+ G L°°(M + ) such that for 
any T > 

lim / sup \U(t,x(t) +y)- U+{t)\ dt = lim / sup \U{t,x(t) + y) - U-(t)\dt = 0. 
n ^°°Jo y e(o,i/n) n ^°°7o s/e(-i/n,o) 

Note that, for each fixed curve, this is equivalent to the convergence for almost every time t. 
Obviously, any BV function verifies this strong trace property. But this requirement is weaker 
than the BV property. Let us emphasize that this notion of trace is more restrictive than the 
strong trace introduced in [43], which is known to be verified for bounded solutions of scalar 
conservation laws. This has been shown in the multidimensional case, first with a non-degeneracy 
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property, in |43j . In the one-dimensional case, a different proof based on compensated compactness 
was proposed by Chen and Rascle [14]. For a general flux function the strong trace problem 
has been solved in the ID case in ,23 . The general multidimensional case has been obtained 
by Panov [351 (33] ( see also Kwon [32], De Lellis, Otto, and Westdickenberg [T7] for interesting- 
generalizations). In the case of systems, this is mainly an open problem. This has been shown 
only for the particular case of isentropic gas dynamics with 7 = 3 for traces in time (traces in 
space can be shown the same way) in [41] . Unfortunately, there are no such results for the strong 
trace property of Definition [1] outside of the usual BV theory. 

Stability of shocks in the class of BV solutions has been investigated by a number of authors. 
In the case of small perturbations in L°° n BV, Bressan, Crasta, and Piccoli [7] developed a 
powerful theory of L 1 stability for entropy solutions obtained by either the Glimm scheme |19] or 
the wave front-tracking method. A simplified approach has been proposed by Bressan, Liu, and 
Yang [5] and Liu and Yang [3JJ. (See also [5].) The theory also works in some cases for small 
perturbations in L°° n BV of large shocks. See, for instance, Lewicka and Trivisa [35J or Bressan 
and Colombo [6j. 

However, our stability result goes beyond the known results valid in the class of BV solutions, 
with perturbations small in BV. Our approach is based on the relative entropy method first 
used by Dafermos and DiPerna to show L 2 stability and uniqueness of Lipschitzian solutions to 
conservation laws [T5J [151 HH] ■ Note that in [T8] , uniqueness of small shocks for strictly hyperbolic 
2x2 systems is shown in a class of admissible weak solutions with small oscillation in L°°r)BV. The 
analysis in [T8] also implies the uniqueness of shocks for 2 x 2 systems in the Smoller- Johnson class 
|38] . In each case genuine nonlinearity is assumed. The ideas of DiPerna were developed further 
by Chen and Frid in the papers |T0l [IT]. In subsequent work, they established, together with 
Li [12) , the uniqueness of solutions to the Riemann problem in a large class of entropy solutions 
(locally BV without smallness conditions) for the 3x3 Euler system in Lagrangian coordinates. 
They also establish a large-time stability result in this context. See also Chen and Li [13] for an 
extension to the relativistic Euler equations. However, no stability in L 2 for all time is included 
in those results. 

Our approach is based on fairly mild assumptions on the system and the Rankine-Hugoniot 
discontinuity. Basically, we need the discontinuity to be extremal (1-shock or n-shock and well 
separated from the other Hugoniot discontinuities), to be either a contact discontinuity or verify 
the Liu condition. In the case of a Liu shock, the corresponding shock curve should satisfy the Liu 
property everywhere (so that the shock speed varies monotonically along the curve), and we need 
also a property of growth of the strength of the shock along the shock curve, where the strength 
is measured via the entropy. Very little constraint is needed on the other shock families. Lax 
properties are typically enough. But we may even relax it to cases where the system is neither 
genuinely nonlinear nor strictly hyperbolic, and even to cases where the shock curves are not 
well-defined. The theory works fine even for large shocks. Note that the present study is another 
step in the program described in [33]. A first step was achieved for scalar conservation laws in 



We will give a precise description of our hypotheses and main results in the next section. First 
let us mention a few particular cases in which our theory applies. Our first examples include the 
isentropic Euler system and the full Euler system for a polytropic gas. Both systems are treated 
in Eulerian coordinates. The isentropic Euler system is the following. 
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We assume a smooth pressure law P : R + — > R with the following properties 

P'(p) > 0, [pP(p)}" > 0. (2) 

As usual, we consider only entropic solutions of this system, namely, those verifying additionally 
the entropy inequality: 

d t r](p, pu) + d x G(p, pu) < 0, 

with 

r,{p, p U )= { -^ + S(p), G(p, pu) = ^gf + pu S'(p), 

and with S"{p) = p^ 1 P'(p) > 0. Note that we need only a single convex entropy, even if in this 
case there exists an entire family of convex entropies. 



The full Euler system reads 

d t p + d x (pu) = 

d t (pu) + d x (pu 2 + P) = (3) 
^ d t (pE) + d x {puE + uP) = 0, 

where E = ^u 2 + e. The equation of state for a polytropic gas is given by 

P = (7 - l)pe (4) 

where 7 > 1. In that case, we consider the entropy/entropy-flux pair 

•q(p, pu, pE) — (7 — l)plnp — pine, G(p, pu, pE) = (7 — l)pu lnp — pu\a. e, (5) 

• i , , P E (P u ) 2 

where, in conservative variables, we have e = 

p 2p z 

For the Euler systems (JT]) and ((3]), we have the following theorem. 

Theorem 1.1. Consider a shock (Ul, Ur) ~ {{pl, ul), (pr, ur)) with velocity a associated to the 
system (pp-©, (resp. (Ul,Ur) — {{pl,ul,El),{pr,ur,Er)) associated to the system ^-^)- 
For any K > 0, there exists Ck > with the following property. For any < s < 1 we have the 
two following cases: 

• // (Ul, Ur) is a 1-shock, then for any weak entropic solution U — (p, u) £ L°°(M. + x K) of 
(resp. U = (p,u,E) G L°°(IR + x R) of verifying the strong trace property of Definition 
\T\ and such that \\(p,u)\\ loo < K, (resp. \\(p, u,E)\\i,^ < K) and 

,0 

\U {x)-U R \ 2 dx<e, / \U (x)-U L \ 2 dx<e A , 

J — OO 

then there exists a Lipschtiz path x(t) such that for any T > we have both 
/ \U(T,x + x(T)) - Ur\ 2 dx <Ce{l + T), \ \U(T, x + x(T)) - U L \ 2 dx < e 4 . 

JO J -oo 
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If (Ul, Ur) is a n-shock, then for any weak entropic solution U G L°°(R + x M) of fl]) (resp. 
of ^Bjj) verifying the strong trace property of Definition^ and such that \\(fl, u)\\l°° < K 
(resp. || (p, u, E)\\l<x> <K), and 



o 
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\U (x) - U R \ 2 dx < e\ / \U (x) - U L \ 2 dx < e, 

J — oo 

then there exists a Lipschtiz path x{t) such that for any T > we have 
/ \U(T,x + x(T)-U R \ 2 dx<e A , / \U(T,x + x(T)) - U L \ 2 dx < Ce(l + T). 



In both case we have 

\x(t)-at\ < Cy/et(l + t). 

Note that this implies both the uniqueness and the stability of the shocks. The solutions are 
not assumed to be away from vacuum. This gives a L 2 stability results up to the translation x(t). 
It is worth noting that the treatment of the entropy on the left and the right of the shock is not the 
same. For example, in the case of a 1-shock, we will show that the total relative entropy on the left 
is strictly decreasing. Hence, if at t = 0, U°(x) — Ul for x < 0, then it stays that way on the left 
of x(t) (while keeping the relative entropy on the right under control). The idea here comes from 
[42] , where a similar stability appeared in the study of a semi-discrete shock for an isentropic gas 
with 7 — 3. For 2x2 systems, all Rankine-Hugoniot discontinuities are extremal, hence Theorem 
11.11 applies to all admissible shocks of (TTJ) . (The theorem also holds for contact discontinuities 
in this case, although the pressure laws are usually nonphysical.) For the full Euler system, all 
shocks are 1-shocks or n-shocks (3-shocks in this case), so again the result applies to any entropy 
admissible shock. However, our result does not provide the stability of contact discontinuities for 
this system. The problem is that the contact discontinuities correspond to 2-waves (with a middle 
eigenvalue). Note that in the isentropic case with P(p) = p 1 (7 > 1), it is enough to assume that 
the initial values are bounded since solutions can be constructed conserving this property (see 
Chen 9J, or Lions Perthame Tadmor [55], for instance). 

We now show an application of our method in the general setting of strictly hyperbolic conser- 
vation laws with either linearly degenerate or genuinely nonlinear characteristic fields. We consider 
an 71 x n system of conservation laws 

d t U + d x A(U) =0, (6) 

which has a strictly convex entropy r\. Assume that A and 77 are of class C 2 on an open state 
domain C K™. We have the following result. 

Theorem 1.2. Assume that the smallest (resp. largest) eigenvalue of V A(V) is simple for all 
V € H, and that the corresponding 1- characteristic family (resp. n- characteristic family) of 0) 
is either genuinely nonlinear or linearly degenerate. Then, for any Vq € £1, there exists K > 
and C > such that, for any entropy- admissible 1-shock or 1-contact discontinuity (resp. n- 
shock or n-contact discontinuity) with speed a and endstates (Ul, Ur) verifying Ul € Bk(Vq) and 
Ur € Bk (Vq), the following is true. For e > small enough, and for any weak entropic solution 
U bounded in Bk(Vq) on (0,T) such that 

/ \U (x)-U R \ 2 dx<e, / \U {x)-U L \ 2 dx<e\ 
J J -00 
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(resp. J °° \Uq(x) — Ur\ 2 dx < e 4 and \Uo(x) — Ul\ 2 dx < e), there exists a Lipschitzian curve 
x(t) such that, for any < t < T , we have both 

/o 
\U(t,x + x(t)) - U L \ 2 dx < e 4 , 

(resp. J °° \U{t,x + x{t)) - U R \ 2 dx < e 4 and j"^ \U(t, x + x(t)) - U L \ 2 dx < Ce(l + t)). In both 

cases we have 

\x{t) - at\ < Cy/et(l + t). 

In particular, this provides L 2 stability, up to a drift, for suitably weak shocks and contact 
discontinuities in a class of perturbations without BV conditions. Note that the assumption 
of genuinely nonlinearity or linear degeneracy applies only to the wave family associated to the 
extremal eigenvalue. No such assumptions are needed on the other wave families. 

The theorems above highlight only a few applications of our theory. In the next section, we 
develop our methods in a more general framework. The assumptions on the Hugoniot curves 
are quite natural and we require no smallness condition on the discontinuities at play. We can 
even relax the strict hyperbolicity condition and consider cases in which the middle eigenvalues 
degenerate and possibly cross each other. 

To conclude this introduction, let us mention that the relative entropy method is also an 
important tool in the study of asymptotic limits to conservation laws. Applications of the relative 
entropy method in this context began with the work of Yau |45) and have been studied by many 
others. For incompressible limits, see Bardos, Golse, Levermore [UI2], Lions and Masmoudi [27] . 
Saint Raymond et al. [201 E3 [32} [36] . For the compressible limit, see Tzavaras [40] in the context 
of relaxation and [H 131133] in the context of hydrodynamical limits. Up to now, this method works 
as long as the limit solution is Lipschitz. It would be of significant interest to extend the method 
to shocks (see [44]). 

2 Presentation of the results 
2.1 General framework 

We want to study a system of m equations of the form 

d t U + d x A(U)=0, (7) 
where the flux function A is defined on an open, bounded, convex set V C M m . 

A:Vci m — ► W m . 

We assume that A e C 2 (V). Additionally, we assume the existence of a strictly convex entropy 

rj : V C K m — > K, 
of class C 2 , and an associated entropy flux 

G : V C E m — ► K, 
of class C 2 , such that the following compatibility relation holds on V. 

m 

djG = ^ dtrjdjAi for any 1 < j < m. (8) 
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If we want to apply our theory to the systems of gas dynamics, we have to define these functions 
on a suitable subset of the boundary of V, namely the points corresponding to vacuum states. For 
this reason, we introduce 

U = V, 

the closure of V. We assume, for simplicity, that A, r/, and G are continuous on U (but with no 
additional regularity up to the boundary) . We denote by U° the subset of U where at least one of 
the functions 77, A, G is not C 1 (typically the vacuum states). 

Remark. It is possible to consider a more general situation in which 77 is unbounded on V. In 
that case, we can add the "vacuum points" in the following way, as in [44] . 

U = {V £ R m I 3V k £ V, lim V k = V, limsup r/{V k ) < 00}. 

We consider here a slightly less general framework, which is enough to treat the Euler systems on 
a state domain with bounded velocities. 

Next, we define, for any V £ V, U £ U, the relative entropy function 
r)(U I V) = r]{U) - r)(V) - Vry(7) ■ (U - V). 
Since 77 is continuous (in fact, convex) on U and strictly convex in V, we have 

ri(U\V) > 0, U eU,V e V, 

and 

rj(U I V) = if and only if U = V. 

Indeed, the following lemma shows that the relative entropy is comparable to the square of the 
L 2 norm. 

Lemma 1. For any compact set Sic V, there exist C\,C<x > such that 

d\U - V\ 2 < r,{U I V) < C 2 \U - V\ 2 , 

for any U € U and V € Q. 

We give a proof of this well-known estimate in the appendix. 

For a pair of states Ul ^ Ur in V, we say that {Ul,Ur) is an entropic Rankine-Hugoniot 
discontinuity if there exists a £ K such that 

A(U R ) - A(U L ) - o-(U R - U L ), 

(9) 

G(U R ) - G(U L ) < a( V (U R ) - V (U L )). 
Equivalently, this means that the discontinuous function U defined by 

{Ul, if x < at, 
Ur, if x > at, 

is a weak solution to ([7]) verifying also, in the sense of distributions, the entropy inequality 

d t r)(U) + d x G{U) < 0. (10) 

Finally, we denote by A~(J7) and \ + {U) the smallest and largest eigenvalues, respectively, 
of VA(U). Hereafter, we assume that A ± (C/) are simple eigenvalues for all U £ V and that 
U -> A ± (J7) lie in L°°(U). (They may be undefined on U°.) 
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2.2 Hypotheses on the system 

First, we assume that for any {U_,U + ) entropic Rankine-Hugoniot discontinuity with ?7_ ^ U + 
we have both ?7_ ^ W° and U + ^14°. (Typically, there is no shock connecting the vacuum.) 

We will consider two sets of assumptions. One set will imply the result on the 1-shock (or 1- 
contact discontinuity), the second set (dual from the first one) will imply the result on the n-shock 
(or n-contact discontinuity). A system satisfying both set of hypotheses, verifies both results. 

First set of hypotheses 

The first set of hypotheses, related to some Ul € V, is the following ((HI) to (H3)). 

(HI) (Family of 1-contact discontinuities or 1-shocks verifying the Liu condition) 

There exists a neighborhood B C V of Ul such that for any U £ B, there is a one parameter 
family of states Su{s) € U defined on an interval [0, su], such that Su(0) — U, and 

A(S u (s))-A(U)=a u (s)(S u (s)-U), s€ [0, 8u ], 

(which means that (U, Su(s)) is a Rankine-Hugoniot discontinuity with velocity ajj(s)). We 
assume that U — > su is Lipschitz on B and both (s, U) — > Su(s) and (s,U) — > <7t/(s) are 
C 1 on {(s, U)\U € 5,0 < s < su}- We assume also that the following properties hold for 
U e B. 

(a) cr[/(s) < for < s < sjj (the speed of the shock decreases with s), and au(0) = \~(U). 

(b) (1-shock) If a'jj ^ 0, then ^-r](U\Su(s)) > (the shock "strengthens" with s) for all s. 

(H2) If (U, V) is an entropic Rankine-Hugoniot discontinuity with velocity a such that V £ B, 
then a > X~(V). 

(H3) If (U, V) is an entropic Rankine-Hugoniot discontinuity with velocity a such that U € B and 
a < \~(U), then (U, V) is a 1-shock. In particular, V = Su(s) for some < s < sjj. 

Second set of hypotheses 

The second set of hypotheses, related to some Ur e V, is the following ((H'l) to (H'3)). 

(H'l) (Family of n-contact discontinuities or n-shocks verifying the Liu condition) 

There exists a neighborhood B C V of Ur such that for every U € B there is a one parameter 
family of states Su(s) G U defined on an interval [0, su], such that Su(0) = U, and 

A(S u (s))-A(U)=a u (s)(S u (s)-U), S e[0, % ], 

(which means that (Sjj(s), U) is a Rankine-Hugoniot discontinuity with velocity crjj(s)). We 
assume that U — > su is Lipschitz and both (s, U) — > Su{s) and (s, U) -4- <tu(s) are C 1 on 
{(s, U)\U G -B, s e [0, s;/]}- We assume also that the following properties hold for U € -B. 

(a) cr^(s) > for < s < sjj (the speed of the shock increases with s), and cr;y(0) = \ + (U). 

(b) (n-shock) If a' v ^ 0, then ^-r\{U\Sxj(s)) > (the shock "strengthens" with s) for all s. 

(H'2) If (V, U) is an entropic Rankine-Hugoniot discontinuity with velocity a such that V G B, 
then cr < X+(V). 
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(H'3) If (V, U) is an entropic Rankine-Hugoniot discontinuity with velocity a such that U £ B and 
a > \ + (U), then (V, U) is an n-shock. In particular, V — Sjj(s) for some < s < su- 

Remarks 

• Note that a given system ([7]) verifies Properties (HI) to (H3) relative to U <E V if and only 
if the system 

d t U-d x A(U)=0, (11) 

verifies Properties (H'l) to (H'3) relative to the same U € V. The properties are, in this 
way, dual. 

• In the case Oy-(s) = for all s, Property (HI) just assumes the existence of a 1-family of 
contact discontinuities. 

• In the case where (Ju(s) is not constant in s, Property (HI) assumes the existence of a family 
a 1-shocks (U, Su{s)) verifying the Liu entropy condition for all s (Property (a)). The only 
additional requirement is (b), which is a condition on the growth of the shock along Su(s), 
where the growth is measured through the pseudo-metric induced from the entropy. This 
condition arises naturally in the study of admissibility criteria for systems of conservation 
laws. In particular, it ensures that Liu admissible shocks are entropic even for moderate to 
strong shocks. Indeed, this fact follows from the important formula 

G(S Ul (a)) - G(U L ) = a Uh (s) ( V (S Ul (*)) - r)(U L )) + f o' Ul (r) V (U L \ S Ul (t)) dr, 

Jo 

which is proved in Section [31 (See also [Tli)l [Ml [3D] .) 

• Hypothesis (H2) is fulfilled under the very general assumption that all the entropic Rankine- 
Hugoniot discontinuities verify the Lax entropy conditions, that is 

k(U-)>a>\ l (U+), 

for any i-shocks (£/_,[/+) with velocity a and any 1 < i < n. Indeed, we need only the 
second inequality, and the fact that Aj(£/+) > \~(U+). 

• Hypothesis (H3) is a requirement that the family of 1-discontinuities is well-separated from 
the other Rankine-Hugoniot discontinuities and do not interfere with them. In the case of 
strictly hyperbolic systems, it is, for instance, a consequence of the extended Lax admissi- 
bility condition 

A l+ i (U+ )> o >A,_ !([/_), 

for all i-shocks (U-, U+), i > 1, with velocity a. Indeed, we use only the second inequality 
and the fact that Aj_i(i7_) > A~ (£/_). Note that we need to separate only the 1-shocks 
issued from B, that is close to Ul- 

• The existence of an entropy 77 implies that the system ([7]) is hyperbolic. Since A £ C 2 (V), 
the eigenvalues of \7A(U) vary continuously on V. In particular, since A ± (J7) are simple for 
U G V, the implicit function theorem ensures that the maps U —> A ± (C7) are in C 1 (V). Note, 
however, that those maps may be discontinuous on U. 
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2.3 Statement of the result 



Our main result is the following. 

Theorem 2.1. Consider a system of conservation laws Q), such that A is C 2 on an open, 
bounded, convex subset V o/R m . We assume that there exists a C 2 strictly convex entropy r] on 
V verifying We assume that rj, A and G are continuous on U = V. Let Ul G V. Assume 
that the system verifies the Properties (Hl)-(HS). Consider U R , G V such that (Ul,U r ) is 
a 1-shock (or 1-contact discontinuity) with velocity a. (This means that there exists s > such 
that Ur = Su L (s) and a = <jxj L {s)). Then there exist constants C > 0, So > such that for any 
< e < Sq and any weak entropic solution U of with values in Lt on (0, T) verifying the strong 
trace property of Definition^ and 

/ \U (x)-U R \ 2 dx<e, / \U (x)-U L \ 2 dx<e\ 

JO J -co 

there exists a Lipschitzian map x(t) such that for any < t < T: 

poo />0 

/ \U{t,x + x{t)) - U R \ 2 dx < Ce(l + t), / \U(t,x + x(t)) - U L \ 2 dx < e 4 . 

JO J -co 

Moreover, 



\x(t) - at\ < Cy/et(l + t). 

This shows that the theory is essentially an L 2 theory. The correction of the position of the 
approximated shock x(t) is fundamental, since the result is trivially wrong without it, even for 
Burgers' equation in the scalar case (see [25]). Part of the difficulty of the proof is to find this 
correct position. Indeed, we will construct it so that /___ T)(U (t, x + x(t))\UL) dx is decreasing and 
remains smaller than e for all time. 

Applying the main theorem on the U(t,x) = U(t, —x), which is an entropic solution to 
we obtain the following corollary. 

Corollary 1. Consider a system of conservation laws such that A is C 2 on a open, bounded, 
convex subset V of R m . We assume that there exists a C 2 strictly convex entropy n onV verifying 
(03). We assume that n, A and G are continuous on 11 — V. Let Ur g V. Assume that 
verifies the properties (H' 1)-(H' 3). Let Ul G V such that (Ul,Ur) is a n-shock (or n-contact 
discontinuity) with velocity a. (This means that there exists s > such that Ul — S\j R (s) and 
a = o~u R (s))- Then there exist constants C > 0, sq > such that for any < e < Eo and any 
weak entropic solution U of ([?[) with values in IA on (0, T) verifying the strong trace property of 
Definition]]^ and 

pco pO 

\Uo{x)-U R \ 2 dx<e\ / \U Q (x)-U L \ 2 dx<e, 

J —CO 

there exists a Lipschitzian map x(t) such that for any < t < T : 

r° 

\U(t,x + x(t))~U R \ 2 dx<e 4 , / \U(t,x + x(t)) - U L \ 2 dx < Cs(l +t). 

J -co 

Moreover, 



\x(t) - at\ < Cy/et(l + t). 
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In particular, if we consider a 2 x 2 system which verifies both (H1)-(H3) for any Ul G V and 
(H'l)-(H'3) for any Ur £ V, then all shocks (and contact discontinuities) are unique and stable 
in L 2 . 

Note that the assumptions on the system are quite minimal. There is an assumption only 
on the wave coming from Ul (or going to Ur in the case of the corollary). There are absolutely 
no assumptions on the other waves (which may not even exist or may be neither genuinely non- 
linear nor linearly degenerate). The extremal property that the shock curve under consideration 
corresponds to the smallest eigenvalue of VA(U) (resp. the largest eigenvalue of VA(U)) is only 
prescribed on a small neighborhood of Ul (resp. of Ur). Finally the theory allows us (via the 
extended set U) to consider weak solutions which may take values U corresponding to points of 
non-differentiability of A and 77. This includes, for example, the vacuum states in fluid mechanics. 

2.4 Main ideas of the proof 

The main idea of the proof is to find a path t — > x(t) for which the total relative entropy of U(t, x) 
with respect to Ul on the left of x(t), given by 



r)(U(t,x + x(t))\U L )dx, 

-00 

is strictly decreasing. The following estimate underlies most of our analysis. 

Lemma 2. If V G V and U is any weak entropic solution of then t](U \ V) is a solution in 
the sense of distributions to 

d t r,{U I V) + d x F(U, V) < 0, 

where 

F(U, V) = G(U) - G(V) - V V (V) ■ (A(U) - A(V)). 

The proof of this lemma is direct from the definition of the relative entropy (Note that V is 
constant with respect to t and x) . From this lemma, and the strong trace property of Definition 
[U we will show that 

- / r]{U(t,x + x(t))\UL)dx<x'(t)T](U{t,x(t)-)\UL)-F(U{t,x(t)-),UL), 
dt J-oo 

for almost every t. Let us fix a small drift velocity vq > of order e. Ideally, we would like to 
have 

F(U(t,x(t)-),U L ) 

x ^=w^mm vo - (12) 

This would preserve the estimate 



rj(U(t, x + x(t)) I U L ) dx < e 4 , 

-00 

which is assumed at the initial time, and ensure that 

, r)(U(t,x + x(t))\U L )dx< -v £ 2 , whenever rj(U(t,x(t)-) | Ul) > e 2 . 
dt J-oo 

If we denote by I the set of time for which rj(U (t, x(t)—) \ Ul) > e 2 , then we must have 

e 4 
\T < " 



v e 2 
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Before going on, let us note that U has little regularity, and (TH21) cannot be solved in the 
classical sense. Hence we can define x(t) only in the Filippov way. (This is roughly the idea. For 
technical reasons, the construction is even more complicated than that.) We will have to check 
carefully that we can do it using only the strong trace property of Definition [T] Even so, we cannot 
ensure that (IT^I) holds almost everywhere. However, we will use the fact that for almost every 
time t, especially when U(t, x(t)+) ^ U(t,x(t)—), the following Rankine-Hugoniot relation holds: 

A(U{t, x(t)+) - A(U{t, x(t)-) = x'{t)(U(t, x(t)+) - U(t, x(t)-)). 

We will make these claims more precise in the sections to follow. 

At this point we need to check that the total relative entropy of U (t, x) with respect to U R on 
the right of x(t) remains under control. A similar estimate to the one above gives 

poo roc 

/ V (U(t,x + x(t))\U R )dx= / r 1 {U Q {x)\U R )dx 
Jo Jo 

+ f {F(U(t,x(t)+),U R ) -x'(t)ri(U(t,x(t)+)\U R )} dt. 
Jo 

Since U is bounded, we have the easy estimate 

J {F(U(t,x(t)+),U R ) - x'(t)r)(U(t,x(t)+) | U R )} dt < Ce. 

In the remaining case, we have to control F(U(t,x(t)+), U R ) — x'(t)r](U(t, x(t)+) \ U R ) for t such 
that r)(U(t,x(t)-) \ U L ) < e 2 , that is 

\U(t,x(t)-)-U L \~e. 
The extremality property of (Ul, U R ) ensures that the worst scenario corresponds to 
U(t,x(t)+) = Su(t,x(i)-)(s), x'(t) = <7p(t, a! (t)_)(s) for given s > 0. 
But, Su(t,x(t)-)( s ) ~ Su L ( s ) + an d the key structural lemma will show that 

F(S Ul (s), Ur) - ct Ul (s) V (S Ul (s) I ITr) < 0, for a > 0. 

Hence, for t £ I, 

F(U(t, x{t)+), U R ) - x'(t)r,(U(t, x(t)+) \ U R ) ~ e. 

The rest of the paper is organized as follows. In the next section, we prove the main structural 
lemmas. In the following section we construct the path t — > x(t). The next one is dedicated to 
the proof of the main theorem. In the last section, we prove theorems on the example stated in 
the introduction. 

3 Structural lemmas 

The first lemma of this section gives an explicit formula for the entropy lost at a Rankine-Hugoniot 
discontinuity ([/_, U + ), where U+ = Su_ (s) for some s > 0. The estimate can be traced back to 
the work of Lax [21] . 
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Lemma 3. Assume (U-,U+) G V 2 is an entropic Rankine-Hugoniot discontinuity with velocity 
a; that is, ({/_, U+) verifies (QJ). Then, for any V <ElA 

F(U+,V) - ar)(U+ | V) < F(U-,V) - a V (U- | V), 

where F is defined as in Lemma\^ Furthermore, if U— £ B, as in Hypothesis (HI), and there 
exists s > such that [/+ = Su_ (s) and a = ajj_ (s) (that is, (f E/+) is a 1- discontinuity), then 

F(U+,V) - a V (U + | V) = F(U-,V) - ar)(U- \ V) + [ o' v (t) V (U- | S v _ (t)) dr. 

Jo 

Proof. Since (£/-, [/+) € V 2 is an entropic Rankine-Hugoniot discontinuity with velocity a we have 
-W V (V) ■ (A(U+) - A{U-)) = -crW V (V) ■ (U+ - [/_), 

and 

G(U+)-G(U-)<a{ri{U + )- V (U-)). 

Summing those two estimates gives the first result. 

Assume now that it is a 1-discontinuity. Then, define 

F^s) = F(Su_(s),V) - F(U^,V), 

F 2 (s) = a(s)(r,(Su_ (s) \ V) - r,{U- \ V j) + f a' v _ (r)r,{U^ \ S v _ (r)) dr. 

Jo 

We want to show that Fx(s) = ^(s) for all s. Since <S[/_(0) = {/_, the equality is true for s = 0. 
Next we have 

F[(s) = ±G(Su_(s))-Vr,(V) ■ ^-A(Su_(s)) 
as as 

= [VtXSW*)) - Vv(V)} ■ ±[A(S V _{8)) - A(U L )}, 

and 

F' 2 {s) = o'u_ (s)[Vrj(V) ■ ((Su_(s) - V) - - V)) - Vn(Su_(s)) ■ (S V _ - U-)] 

+ o v _ {a)\yr,(Su_ (s j) - V V (V)} ■ S'^ (s) 

= [VviSuAs)) - V V (V)} ■ £[<ru-(a)(Su-(*) - U L )\. 

as 

Using the fact that (t/— , Su_ (s)) with velocity ajj_ (s) verifies the Rankine-Hugoniot conditions, 
we get 

F[(s) = F' 2 {s) fors>0. 

□ 

The next lemma is fundamental to the theory. It is a slight variation on a crucial lemma of 
DiPerna p] . 

Lemma 4. For any U £ B and any s, so € [0, su], we have 

F(Su(s), Su(s )) - o-u(s)r){Su(s) \ Sv(s a )) = f o f v (T){n(JJ\S v (T)) - V (U \ SuM)) dr < 0. 

J so 

Proof. We use the estimate of Lemma [3] twice with V = Sjj(so) and U- — U . The first time we 
take U+ = Sjj(s), and the second time U+ — Su(sq). The difference of the two results gives the 
lemma. □ 
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4 Construction of t — > x(t) 



In this section, we contruct the path t — > x(t) and study its properties. Unlike the construction 
in [25] . we do not appeal directly to the theory of Filippov. Instead, we build x(t) as the limit of 
approximate solutions, and then show that the Filippov properties are satisfied. The procedure is 
interesting in its own right and has similarities with the method of Hermes (see [21] )• We will focus 
on the family of 1-discontinuities. The corresponding results are easily obtained for the family of 
n-discontinuities. 

Throughout this section, we assume that (Ul, Ur) € V 2 is a fixed 1-discontinuity with velocity 
<t, and that U is a fixed weak entropic solution of ([7]) verifying the strong trace property of 
Definition [TJ We also fix e > 0. As before, X~(U) denotes the smallest eigenvalue of VA(U) and 
is defined for all U £ U \ 14°. We define on U the velocity function 



V(U) 



\f{u,u l ) E>x -^_ E i if u eu\(U°u{u L }), 



U(u\u L ) 

\-(U L )-e, if U = U L , (13) 

F (U,U L ) _ ; C rr^,0 



I v(u | U L ) 
We have the following lemma. 



e, if UGU°. 



Lemma 5. The function U — > V(U) is continuous on 14 \ IA° , and is bounded and upper semi- 
continuous onU. Moreover, there exists a constant C (independent of e) such that if U € B and 
r/(U\U L ) < e 2 , then 

V(U) > \-{U L )-Ce. 
Note that this map may not be continuous on U°. 

Proof. First note that the function U — > ^jt^i is continuous onW\ {Ul}- Also, by assumption, 
A~([/) is continuous onU\U Q and bounded onW. So U — > V(U) is continuous onU\ ({(7l}UW ), 
and upper semi-continuous onW\{(/i}. At the point Ul, we have 

VuF(U,U l )\ u=Ul =0, VIF(U,U l )\ u=Ul = D 2 n(U L )WA(U L ), 
Vuti(U\U l )\ u=Ul =Q, ^t,{U\U l )\ u=Ul =D 2 t,(U l ). 

Owing to the strict convexity of r\ at Ul , an expansion near U l gives 

F(U, U L ) (U - U L ) T D 2 n(U L )WA(U L )(U - U L ) 



r,(U | U L ) " (U - U L ) T D 2 V (U L )(U - U L ) 



0{\U-U L \). (14) 



Since D 2 i](Ul) is symmetric positive definite and D 2 t](Ul)VA(Ul) is symmetric, those two ma- 
trices are diagonalizable in the same basis. This gives 



Therefore, 



which implies 



\-{Ul)D\{Ul) < D 2 n{U L )VA(y L ) < A+ (U L )D 2 i 1 {U L ). (15) 



hminf F ^ T f TJ L \ -e>\-(U L )-e= V(U L 



lim VOJ) = \-{U L )-e. 

U^U L 
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So, V is continuous at U = Ul and bounded on IA. 
More precisely, (fl"4l and (fT5)) give that 

Hffj>A-(%)-W-^l). 
Since [/ — > A~([/) lies in C 1 (B), we have also that 

x-(u)>\-(u L )-o(\u-u L \). 

This provides the last statement of the lemma. □ 
For any Lipschitzian path t — > x(t) we define 

W*) = max{V(t/(*,z(i)-)), *(*)+))}, 



V min (t) = 



min {F(f7(t, ar(t)— )), V(U(t, x(t)+))} , if both U(t, x(t)-), U(t, x(t)+) <£ U°, 
— oo, otherwise. 



This section is dedicated to the following proposition. 

Proposition 1. for any (Ul, Ur) G V 2 1- discontinuity with velocity a, e > 0, and U a weak en- 
tropic solution of £?]) verifying the strong trace property of Definition^ there exists a Lipschitzian 
path t — > x(t) such that for almost every t > 

V miD (t) < x'(t) < V max (t). 

Proof. Consider the function 



S,x)= / V(U(t,x+^)) 
Jo n 



dy. 



By virtue of Lemma [5l v n is bounded, measurable in t, and Lipschitz in x. We denote by x n the 
unique solution to 

f x' n (t) = v n (t,x n (tj), t > 0, 
\ x n (0)=0, 

in the sense of Caratheodory. Since v n is uniformly bounded, x n is uniformly Lipschitzian (in 
time) with respect to n. Hence, there exists a Lipschitzian path t — > x(t) such that (up to a 
subsequence) x n converges to x in C°(0,T) for every T > 0. We construct V max (i) and Knin(*) as 
above for this particular fixed path t — > x(t). Let us show that for almost every t > 

lim [x' n (t)-V max (t)] + =0, (16) 

n—too 

lim [Knn(«) -a£(t)]+ = 0. (17) 
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Both limits can be proved the same way. Let us focus on the first one. We have 



K(t) - v max (t)}+ 



V(U(t,x n (t) + ^-))dy-V max {t) 
n 



[V(U(t,x n (t) + ?-)) -V max (t)]dy 



< 



V(U(t,x n (t) + ^))-V m ^{t) 



dy 



< ess sup 

1/6(0,1) 



V(U(t,x n (t) + ^))~V m ^(t) 



< esssup [V(U(t,x{t)+z))-V max (t)] 

where, for a given t > 0, e n —> is chosen so that x n (t) — x(t) G (— £„,£„ — — ). We claim that 
for almost every t > 0, the last term above goes to zero as n — > oo. Indeed, fix t > for which 
U(t,x(t) + •) has a left and right trace in the sense of Definition [TJ That is, 

lim \ esssup \U(t,x(t) + y) - U + (t)\ } = lim \ ess sup \U(t, x(t) - y) - U-(t)\ } = 0, 



e-s-0 



3/e(0,e) 



e->0 



Since V is upper semi-continuous on we have that for all r > there exists 6 > such that 
\U-U±(t)\<6 => [V(U)-V(U±{t))]+<r. 

Therefore, 

lim \ esssup [V(?7(t,a;(t) ±j/)) - V(C/± (*))], i =0, 



e->0 



ye(o,e) 



and it follows easily that 



lim <^ esssup [V(U(t,x(t) + z)) - V max {t)}_ 



2£(-E,e) 



This verifies the claim above and finishes the proof of (fT6|). The proof of (|17p is similar; we just 
use the continuity of V on U \ U° and the definition of V m - ln for U Gli . 

Finally, the sequence x' n converges to x' in the sense of distributions. Also, the function [•]+ is 
convex. Therefore, integrating (fT6|) and (|T7| on [0, T] and passing to the limit, we obtain 

[V min (t) - x'(t)}+ dt= [ [x'(t) - V max (t)}+ dt = 0. 



'0 JO 

In particular, for almost every t > we have 

V min (t) < x'(t) < V m3X (t). 



□ 
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We end this section with an elegant formulation of the Rankine-Hugoniot condition and related 
entropy estimates, as originally presented by Dafermos in the BV case. We show that the estimates 
remain true for solutions having the strong trace property (in fact, the strong trace property defined 
in 43 suffices) . The proof is given in the appendix. 

Lemma 6. Consider t — > x(t) a Lipschitzian path, and U an entropic weak solution to ^ verifying 
the strong trace property. Then, for almost every t > we have 

A(U(t, x(t)+)) - A(U(t, «(*)—)) = x'{t)(U(t, x(t)+) - U(t, x{t)-)), 
G(U(t,x(t)+)) - G(U(t,x(t)-)) < x'(t)(r,(U(t,x(t)+)) - V (U(t,x(t)-))). 

Moreover, for almost every t > and V G V 

j f J v(U(t, y + x(t)) \V)dy< -F(U(t, x(t)-),V) + x'(t) V (U(t, x(t)-) \ V), 

d r°° 

- / r](U(t, y + x(t)) \V)dy< F(U(t, x(t)+),V) - x'(t)n(U(t, x(t)+) \ V). 



dt J 

5 Proof of Theorem 12.11 

This section is dedicated to the proof of our main result, Theorem 12.11 We consider a system of 
conservation laws (JTJ) such that A is C 2 on a open, bounded, convex subset V of R m . We assume 
that there exists a strictly convex entropy 77 of class C 2 on V verifying ©. We assume that rj, 
A and G are continuous on U = V. Let (Ul,Ur) G V 2 be an entropic 1-discontinuity. For the 
neighborhood B 3 Ul given by (HI) we define 

4 = \^v{u\u L ). 

Consider U weak entropic solution of with values in U on (0, T) verifying the strong trace 
property of Definition [TJ and 

/o 
v(U (x)\U L )dx<e 4 , 
-00 

for some e < £o- Consider the path t — > x(t) constructed in Proposition!]] We first prove a pair of 
lemmas which provide estimates for the rate of change of the total relative entropy with respect 
to x{t). 

Lemma 7. For almost every time < t < T , we have 

jj v(U(t,y + x(t))\U L )dy<0. 

Furthermore, if we let X C (0, T) denote the set of time such that n(U(t,y + x(t)—) | Ul) > £ 2 , 
then we have 

|Z| < E. 

Proof. From Lemma [SJ we have for almost every t > 

j t f° V (U(t, y + x(t)) I U L ) dy < -F(U(t, x(t)-), U L ) + x'(t)n(U(t, x(t)-) \ U L ). 



1G 



At times t for which U(t, x(t)—) = Ul, the right-hand side vanishes (and obviously t ^ I). For 
the other times t, Proposition [T] ensures that either 

A*) < v(u(t, *(*)-)) < ^^Ey^ ~ £ < or At) < v(u(t, x(t)+)). 

For the times corresponding to the first case, we have 

j t J v(U(t, y + x(t)) | U L ) dy < -e V (U(t, x(t)-) \ U L ) < 0. 
Especially, for every such time in X, we have 

jj ri(U(t,y + x(t))\U L )dy<-e 3 . (18) 

For the second case, we appeal to Lemma 6. For almost every time we have either U(t,x(t)—) = 
U(t,x(t)+) (which is covered by the first case) or x'(t) is the velocity associated to an entropic 
Rankine-Hugoniot discontinuity (U(t,x(t)—),U(t,x(t)+)). If U(t, x(t)+) e B C V, we have 

x'(t) < V(U(t,x(t)+)) < \- (U(t, x(t)+)) -e, 

which is in contradiction with Hypothesis (H2). Hence, for almost every time t not covered by the 
first case, U(t,x(t)+) ^ B, and so r)(U(t,x(t)+) \ Ul)> £q- In that case, we deduce from Lemmas 
[3] and ©that 

j t J° V (U(t, y + x(t)) | U L ) dy < -F(U(t, x(t)-), U L ) + x'(t) V (U(t, x(t)-) \ U L ) 

< -F(U(t, x(t)+), U L ) + x'(t)r)(U(t, x(t)+) | U L ) 
<-e V (U(t,x(t)+)\U L ) <-e 3 . 

In particular, for every t S I, we have (|18p . Integrating in time, we find 



\Z\e 3 < / v(U Q (x)\U L )dx<e 4 

J — oo 



which gives the desired result. □ 

Before we prove the second lemma, we mention the following basic fact. 

Remark. If / : U — > M is C 1 on V and bounded on U, then for any fixed V E V, there exists 
C v > such that for all U € U 

\f(U)-f(V)\<C v \U-V\. 
The proof of this fact follows the same outline as the proof of Lemma [T] 
Lemma 8. There exists a constant C > such that for every < e < Eq 
F(U(t, x{t)+), U R ) - x'(t)n(U(t, x(t)+) | U R ) < Ce 

for almost every t ^ I. 
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Proof. First note that t $ X implies U(t,x(t)-) e B' = {U 6 B \ <q(U\U L ) < el) C B. Also, 
from Lemma [1] we have 

\U(t,x(t)-)-U L \ < Ce. 

We consider three cases. 

1. (Points of continuity.) For almost every t ^ 1 such that 

U(t,x(t)+) = U(t,x(t)-), 
we have, from Proposition [T] and Lemma [5l 

x'(t) > V(U(t,x(t)+) > X~(U L ) - Ce. 
Hence, changing the constant C from line to line if necessary, 

F(U(t, x(t)+), U R ) - x'(t)r,(U(t, x(t)+) | U R ) 

< F(U(t, x(t)-), U R ) - \-(U L WU(t, x(t)-) | U R ) + Ce 

< F(U L , U R ) - \-(U L )r](U L I U R ) + Ce 

< F(S Ul (0), S Ul (s )) - <j Ul (0)v(S Ul (0) | S Ul (s )) + Ce < Ce, 
where, in the last inequality, we used Lemma 2] with Su L (so) = Ur. 

2. (Discontinuities excluding 1-shocks.) If (U(t,x(t)—),U(t,x(t)+)) is an entropic Rankine- 
Hugoniot discontinuity with velocity x'(t) > A - (U(t, x(t)— )), then 

F(U(t, x(t)+), Ur) - x'(t) v (U(t, x(t)+) | U R ) 

< F(U(t, x(t)-), Ur) - x'(t) v (U(t, x(t)-) I Ur) 

< F(U(t, x(t)-), Ur) - \-(U_) V (U(t, x(t)-) I Ur) 

< F(U L , U R ) - A- (U L ) V (U L I U R ) + Ce< Ce. 

In the second line we used Lemma [3l and in the fourth line Lemma |4j 

3. (1-shocks.) On the other hand, if (U(t,x(t)—),U(t,x(t)+)) is an entropic discontinuity 
with velocity x'(t) < A~(U(t,x(t)—)), then, by Hypothesis (H3), it is a 1-shock and, since 
U — > su is Lipschitz on B, there exists s' £ [0,su L ] such that we have both U(t,x(t)+) = 
Su(t.x(t)-)(s) = S Ul (s') + 0(e) and x'(t) = o-u(t,x(t)-)( s ) = a u L ( s ') + 0(e). As noted 
above, U(t,x(t)—) e B' C B. Moreover U — > su is continuous on B' , which is closed, and 
there is no state U+ € U° such that (U, U+) is an admissible discontinuity. Hence the set 
{Su(s) | U £ B', s € [0, su]} is at a positive distance from U° (where rj and A may be not 
C 1 ). The functions F(-,Ur) and r](- 1 Ur) are then uniformly Lipschitz on this set. Hence, 
by Lemma |4l 

F(U(t, x(t)+), Ur) - x'(t)r,(U(t, x(t)+) \ Ur) 

< F(S Ul (s'), Ur) - o- Ul (s') V (S Ul (s') I U a ) + O(e) < O(e). 



□ 

We can now finish the proof of the theorem. Lemma [7] implies that for all < t < T 

V (U(t,y + x(t))\U L )dy< f V (U (y)\U L )dy<e 4 . (19) 
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On the other hand, 

r)(U{t,y + x(t))\U R )dy 

r,(U (y)\U R )dy 



(p,t)ra dr 



— { / r)(U (r, y + x(r)) \ U R ) dy \ dr 



'(o,t)ni 

Since F, x', and ij are bounded and \1\ < e, we have 



— |jf 7?(J7(r,y + a;(r))|f7ii)dyJ> dr 



i uo 



rj{U(r, y + x(r)) \ U R ) dyj dr 

\F(U(r,x{r)+),U R )-x / (r)r 1 (U(r,x(r)+) | U R )\ dr < C\l\ < Ce. 
Applying this estimate together with Lemmas [T] and [8] above, we obtain 

V (U(t, y + x(t)) | U R ) dy < Ce(l + t). (20) 

It remains to show that \x(t) — at\ < C\J et(\ + 1). We denote by U± the function defined by 

U±(x) 



U L , ifx<0, 
U R , if x > 0. 



Note that the exact shock solution is U±(x — at). By virtue of (fT9|) and (|20|) . we have for every 
< t < T 

\\U(t, •) - U ± (- - x(t))\\ L 2 < Cy/(l + t)e. 

for some C > "J~~^, where G\ is given by LemmaHJ Let M = \\x' (t)\\ . Then for T > 0, we 
consider an even cutoff function 6 € C°°CR) such that 



4>(x) = l, 
4>{x) = o, 

cf>'(x) < 0, 



if \x\ < MT, 
if \x\ > 2MT, 
if x > 0, 



< 2(MT)" 1 , ifxeM. 
Then, for almost every < t < T, we have 







[d t U + d x A(U)} dxdt 



(x) [U±(x - x(t)) ~ U±(x)} dx - 



A(U±(x — x(s)))4> (x) dx ds 



+ / <p(x)[U{t,x)-U±(x-x(t))] dx+ I <f>(x) [U±(x)-U°(x)] dx 



[A(U(s, xj) - A(U±{x - x(s)))] tj>'{x) dxds. 
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The terms on the second line above to reduce to 

x{t)(U L - U R ) - t(A(U L ) - A{U R )) = (x(t) - at)(U L - U R ). 

The third line can be controlled by 

UWl^V) {\\U(t, •) - U ± (- - x(t))\\ LHm + \\U°- U ± \\ L 2 (m ) < CVIMT^(1 + T)e. 

Finally, since A has a suitable Lipschitz property at the points Ul,Ur € V (see the remark 
preceding Lemma [5]), the last term has the following bound: 

[A(U(s, x)) - A(U ± (x - x(s)))} (j)'(x) dxds 



t r2MT 



< 



2C 



MT J 

Combining the estimates above we obtain 



\U(s,x) — U±(x — x(s))\ dxds 

J -2MT 

V4MT\\U(s, •) - U±(- - x(s))\\ L 2 m dxds < Cy/T(l + T)e. 



W) - at\ < C ^ u T l {1 _ + u ^ £ < C^WVTTs. 



This concludes the proof of the theorem. 



6 Proof of the theorems stated in the introduction 



We can now prove Theorems 11.11 and 1 1 . 21 which were stated the introduction. In each case we want 
to show that the hypotheses of Theorem 12.11 are fulfilled. 



6.1 Isentropic Euler system 

To illustrate the foregoing ideas, we apply our results to the isentropic Euler system (TTJ). We will 
show that our hypotheses are satisfied for a fairly large class of pressure laws P{p). 
First, let us take 

V = {(p,pu) e 1+ x 1 | < ||(/3,it)||ioo < K) . 

The petal-shaped region V is open, bounded, and convex. In particular, u is bounded in V, which 
ensures the continuity of r\, A, and G at the vacuum state (0, 0) G U°. (Here, the set U° contains 
a single point.) One can check that differentiability fails at this point. 
The eigenvalues for this system are 



Xl(p,pu) =U- y/P'{p), \ 2 (p,fm)=U+y/P'(p). 



(21) 



To ensure hyperbolicity we assume P'(p) > 0. If {a, {pl, ul), (pr, ur)} satisfies the Rankine- 
Hugoniot condition, then the following relations hold 



i PR 
V PL 



Pr - Pl 

PR- PL J 



(u L - ur) 2 = 



V PR 



Pr-Pl 



PR - PL 



(Pr-Pl)(pr-Pl) 

PRPL 



(22) 
(23) 
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In view of (|2"Tj) . the plus and minus signs in (|2"2"j) distinguish the 2-shocks (pu < pl) and 1- 
shocks (pa > Pl), respectively. It will not be necessary to distinguish between shocks and contact 
discontinuities here, so we will call all such discontinuities shocks. Let us assume that the entropy 
admissible shocks are characterized by the relation mj, > ur. The opposite case is treated similarly. 
With respect to the hypotheses (HI) and (H'l), we take 



S(p,pu)(s) = P + s, pu±p< 



! (P(p + S )-P(p))s 



p(p + s) 



C(p,pu)(s) =U± 



P(p + s)-P(p) 



where s > represents the change in density, in absolute value, across the shock. The plus and 
minus signs correspond to the cases (H'l) and (HI), respectively. Note that, in both cases, the 
functions (s, U) — > Su{s) and (s,U) —> <jjj{s) are C 1 on K+ x V. Also, it is straightforward to 
verify that the function U —¥ su, where 

s v = sup {s G R + | Su(s) G V for all s G [0, s ]} , 

is Lipschitz on B C V, provided the neighborhood B (corresponding to a fixed left or right 
endstate) is small enough. In particular, the terminal value sjj is finite and does not correspond 
to the vacuum state. 

Using only the hyperbolicity assumption on P, we can verify Hypotheses (Hlb) and (H'lb). 
Indeed, the relative entropy in this case is given by 

P, 



V((P, pu) | (po, PoUq)) = - u ) 2 + S(p\ p ). 



Therefore, 



^((p,HI% pu) (s)) = ^{f 



(P(p + s)-P(p))s 



s{p\p- 



2(p- 



p(p + s) 

[s(p + s)P'(p + s)+ p(P(p + s)- P(p))} + sS"(p + s) > 0. 



Finally, observe that the global Liu-admissibility hypotheses, (Hla) and (H'la), are satisfied 
if and only if the function 

P(p + s)-P(p) 



is nondecreasing. The derivative of 4> can be expressed as follows. 



l 

ps 2 



p+s 



(q-p)[qP(q)}"dq if S ± 0; 



2p 



if 



Therefore, the Liu condition holds if and only if 

\pP(p)\" > o. 
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This may include, for instance, systems which fail to be genuinely nonlinear. On the other hand, 
it is interesting to note that the genuine nonlinearity condition, \pP(p)]" + 1 0, holds if and only 
if a' 7^ 0; that is, if and only if the shock speed is strictly monotone along each shock curve. In 
fact, the equivalence of genuine nonlinearity and the global strict Liu-condition can be shown for 
a wide class of conservation laws, as suggested by Liu [25]. Finally, it follows from (|22l) and the 
analysis above that the Lax admissibility conditions also hold globally when \pP(p)]" > 0. This 
verifies (H2)-(H3) and the corresponding dual conditions. Therefore, Theorem 12 . 1 1 applies and the 
proof is complete. 



6.2 Full Euler system 

Now let us show that our theorem applies to the full Euler system for a polytropic gas. The idea is 
generally the same; however, our presentation is less explicit due to the complexity of the system. 
First, we introduce the state space 

V = {{p,pu,pE) G R + x R x R+ I < \\(p,u,E)\\ L «, < K) . 

As before, the region V is open, bounded, and convex; and the functions rj, A, and G are continuous 
at the vacuum state (0, 0) G U°. 

To simplify the remainder of the presentation, we work with the nonconservative variables 
(p, u, e). The sound speed in the case of a polytropic pressure law (|4]) is given by 



d p P + p~ 2 Pd e P = ^7(7 -l) e . 



Hence the eigenvalues of the system are 



Ai(p, u, e) =u - \Ay(7- l)e, \ 2 (p,u,e)=u 

consider a discontinuity {cr, (pl,ul, eL 
tion. Excluding contact discontinuities (ul 

fiK /7-1 7 + 1 \Pr] 

U L ±\ 4/— h— — = cr = u R ± 



A 3 (p, u, e) = u + y/j(j-l)t 



(24) 



We consider a discontinuity {a, (pl,ul,6l), (pR,,UR,en)} verifying the Rankine-Hugoniot condi- 

ur), the following relations hold. 



PL 



27 



27 



Pl 



PR 



'7 — 1 7 + I 



27 





7 + 


1" 


PR _ 1 
PL 


Pr 


7- 


1 


Pl 


"7- 


- 1 


_ PR 




.7" 


- 1 


PL 



C J L 



Pr Pl_ 
Pl Pr 



7 - 1 PR _ 

7 + 1 PL 7 



27 

7 + 1 



Pl 
Pr 



1 



(25) 



(26) 



We refer the reader to [35] for a proof of these relations. Note that the plus and minus signs in 
l|25j) distinguish the 3-shocks and 1-shocks, respectively. Let us fix a left endstate (pl,ul,gl) £ 
(K + x K x R + ). The entropy admissible shocks can be parametrized by any of the three quantities 
Pr> Pr-, and en, which are related by (|26|). The following conditions are equivalent and characterize 
the family of 1-shocks. 

Pr > Pl, Pr > Pl, e R > e L . (1-shock) 
Similarly, the 3-shocks are characterized by 

Pr<Pl, Pr < PL, e R < e L . (3-shock) 
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Note carefully that, for p^ fixed, pr is bounded above and bounded away from zero, whereas Pr 
and eji range over all positive real numbers. Furthermore, each of the quantities is an increasing 
function of the others as indicated above. While we omit explicit formulas for S^ p _p U p E) (s) and 
a (p,pu,pE)i s ), the regularity properties of (HI) and (H'l) are easily satisfied. 

Note that ([M]) - ([2^|) easily imply that both shock families verify the Lax and Liu admissibility 
conditions globally. Therefore, it remains only to check that the relative entropy increases along 
the shock curves. We will prove that (Hlb) holds for the family of 1-shocks. The proof of (H'lb) 
for the family of 3-shocks is essentially the same. A somewhat tedious calculation shows that 

v((pl,PlUl,PlE l ) I (p R , prur, prE r )) = (7 - l)h(p L | p R ) - p L ln-(e L | e R ) + - u R ) 2 , 

2e R 



where h(x) = xlnx. Observe first that h(x) and —In a; are both strictly convex functions. Since 
IP-R — Pl\ and \cr — &l\ are increasing along the shock curve, the first two terms in the formula 
above are also increasing. To deal with the last term, we recall (|23|) which also holds for the full 
Euler system and which we rewrite as follows. 



(u L - u R ) = 

PR 



1-^ 
Pr 



PL 



Hence we have 

PL 

2e R 



(u L - UrY 



(7 - ^IPlPr 
2Pr 



(u L - u R ) 2 = -(7 - l)p L 



1 



Pl 
Pr 



PR 
PL 



1 



Finally, we note that the bracketed terms are both positive and increasing along the 1-shock curves, 
hence the product is increasing. This establishes hypothesis (Hlb), and the theorem follows. 



6.3 Strictly Hyperbolic system with small amplitude 

The proof of Theorem 11.21 follows the same outline as the proof of Theorem 12.11 We simply need 
to check that our structural hypotheses are verified locally. Again, we restrict our attention to 
1-shocks and 1-contact discontinuities. The argument for n-shocks and n-contact discontinuites is 
similar. 

By assumption, the state space il is open and X~(V) is a simple eigenvalue for all V £ O. 
Therefore, in a neighborhood of each state there exists a 1-shock curve Sy (s), of class C 1 , 

verifying 

Ai(Sy(*)). (27) 

If the 1-characteristic family is genuinely nonlinear, the parameterization of Sv(s) is chosen so 
that the right-hand side of (f27| is negative. For linearly degenerate fields, the right-hand side 
vanishes and <Jy{s) is identically zero for all s along the shock curve Sv(s). In either case, a 
straightforward continuity argument shows that for each Vq E there exists K > (sufficiently 
small) so that for all V £ Bk{Vq), the 1-shock curve Sy{s) exists in a neighborhood of V which 
contains Bk{Vq) and such that <j' v {s) < for Sv(s) £ Bk(Vo). This means that the Liu condition 
(Hla) is valid with respect to any fixed Ul £ Bpc(Vo). In terms of the notation of Section |2~21 we 
have taken V = B = B K (V ). 

Now, with respect to property (Hlb), an easy calculation shows that 

4- V(U I Su(s)) = S'^sf ■ D 2 v (Su( S )) ■ [Su(a) - U]. 
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Since D rj is positive definite, the previous quantity is nonnegative for s > sufficiently small 
(since Sjj(s) — U and S' u (s) T point in nearly the same direction). Hence, (Hlb) holds for all 
Ul £ Bk(Vo), provided K is small enough. 

For genuinely nonlinear fields, it is well-known that the Lax and Liu entropy conditions are 
equivalent for weak shocks. Therefore, it suffices to check conditions (H2) and (H3) among all en- 
tropic Rankinc-Hugoniot discontinuities (£/_,[/+) excluding the family of 1-discontinuities. Keep 
in mind that, in each case, both U- and U+ are contained in the closure of Bk(Vo). Since all 
discontinuities under consideration are weak, we have 

<t(U-,U+) = \{U±) + 0(\U+-U-\), 

where X(U±) is some intermediate eigenvalue of VA(U±). Since WA is continuous and Ai = A~ 
is a simple eigenvalue, there exists 8 > such that X~(V) < X(V) — S for all V £ Bk(Vo). Then, 
for K small enough, it follows that a(U-,U+) > X~(U±) for all (intermediate) Rankine-Hugoniot 
discontinuities with endstates U-, U+ £ Bk{Vq). This completes the proof. 

A Proof of Lemmas [I] and [6] 

We first give the proof of Lemma [TJ 
Proof. On V 2 , we have 

1 p i 

T r>2„ 



r)(U\V)= / (U-V) ■ D 2 r/(V + st(U -V)) ■ (U -V)tdsdt. 
Jo Jo 

Consider a convex, compact set (l such that fi C SI C V and Sl c n fl = 0, where il c — V \ il. Also, 
let < < A^ < oo denote, respectively, the smallest and largest eigenvalues of D 2 r\ on Cl. 

Then, for any U £ and V £ we have 

A~ A+ 
-Y\U- V\ 2 < rj(U | V) < \U - V\ 2 . 

On the other hand, there exists a constant C > such that 

inf rj(U\V)>-, sup r)(U \ V) <C. 

Indeed, rj(- \ •) is continuous on the compact set fl c x J7. So it attains its maximum and minimum. 
The maximum is bounded. The mininum is attained at some point (Uq, Vq) with Uq ^ Vq, and so 
v{Uo I Vq) 7^ 0. Hence the results holds true with 



f i : min { — p-, —, p7 } . ("'■> mux ! 

^ 2 C sup \U ~V\ 2 ' - ' 

uen c ,ven 



We now give the proof of Lemma [6l 



At C 



2 ' inf \U-VY 



□ 
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Proof. Since A, rj, and G are continuous on U, the functions of (t,x) A(U), rj(U \ Ul), v(U \ Ur), 
F(U, Ul), and F(U, Ur) also verify the strong trace property of Definition [TJ 

Consider the family of mollifyer functions 

4>e{y) = -<f>i(~), 
e e 

where <pi G C°°(R) is nonnegative, compactly supported in (—1,0), and such that J 4>i(y) dy = 1. 
We define 

poo 

$e(aO = / 4>e{y)dy. 

J X 

Note that $ e is bounded by 1. It is equal to 1 for x < —e and equal to for x > 0. Especially, it 
converges to l{a;<o} when e goes to 0. 



We consider also a test function in time only ip g C°°(R + ), tfi > 0, compactly supported in 
(0,T) C (0,oo). 

For any function V verifying the strong trace property of Definition [TJ we have 

iP(t)<t> e (y - x(t))V(t,y) dydt 



i>(t)V(t,x(t)-)dt+ / / t/;(t)<l> e (y-x(t))[V(t,y)-V(t,x(t)-)}dydt. 
o Jo 



The last term is smaller than 



IMU- I I Mv)\V(t,x(t)-ey)-V(t,x(t)-)\dydt 
Jo Jo 

<Wl-II^iIU- sup f \V(t,x{t)-y)-V(t,x(t)-)\dt. 

y€(0,e)Jo 

This converges to for e going to 0, thanks to the strong trace property. Hence 



lim/ iP(t)4> e (y-x(t))V(t,y)dydt= i/>(t)V{t,x(t)-)dt. (28) 
e ^°Jo Jr Jo 

In the same way, we show that 

pOO p poo 

lim / i,{t)<t> e {x{t)-y)V{t,y)dydt= x/j(t)V(t,x(t)+) dt. (29) 
Jo Jr Jo 

And so 

pOO p poo 

lim / Hj(t)(<j> e (y-x(t))-ct> e (x(t)-y))V(t,y)dydt= / ii(t)[V(t,x(t)-) - V(t,x(t)+)} dt. 
s ^°Jo Jr Jo 

(30) 

We consider the test function in time and space 

rp(t)$ e (x-x(t)). 
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For the equation of Lemma [5] with V G V, we get: 

ip{t)$' e (x - x(t))[x'(t)r)(U(t, x)\V)- F(U(t, x), V)} dxdt 
>- I / i(>'(t)<f> e (x - x(t))ri(U(t,x)\V)dxdt. 



JO JR 

Passing to the limit and using (f28f , we find 

roo />0 poo 

- ip'(t) r]{U{t,x{t)+x)\V)dxdt < ^(t)[a/(t)T)(U(t,x(t)-)\V)-F(U(t,x(t)-),V)]dt. 

Jo J-oo JO 

This is the desired result in the sense of distribution. 

Wc consider now the test function in time and space 

1>(t)* e (x(t) - x). 

For the equation of Lemma [2] with V G V, we get in the same way, using (|2T))) . 

/>oo />oo 

ip'{t) / r)(U(t,x(t)+x)\V)dxdt< / ip(t)[-x'[t)r](U[t,x(t)+) \ V)+F(U(t,x(t)+), V)} dt. 
o Jo Jo 

We take, now, as a test function 

#)[$e(z - x(t)) + $ e (a;(t) - s) - 1]. 
Note that this function converges to in L 1 when e converges to 0. For equation (JT]), this gives 



o 



ip(t)[& e (x - x(t)) - <f>' e (x(t) - x)][x'(t)U(t,x) - A(U(t,x))]dxdt 
ip'(t)[$ e (x - x{t)) + $ e (x(t) -x)- l]U(t,x)dxdt. 
The right hand side term converges to 0. Thanks to (|30l) . the left hand side term converges to 
<tp(t)[x'{t)(U(t, x(t)+) - U(t, x{t)-)) - (A(U(t, x(t)+)) - A(U(t, x{t)-)))] dt = 0. 



This provides the first equality in the sense of distribution. The second one can be proven the 
same way from (|10[) . □ 
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